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ABSTRACT 

We present Spitzer measurements of the aromatic (also known as PAH) features for 35 Seyfert 
galaxies from the revised Shapley-Ames sample and find that the relative strengths of the features 
differ significantly from those observed in star-forming galaxies. Specifically, the features at 6.2, 
7.7, and 8.6 fxva are suppressed relative to the 11.3 feature in Seyferts. Furthermore, we find 
an anti-correlation between the L(7.7 //m)/L(11.3 /im) ratio and the strength of the rotational H2 
emission, which traces shocked gas. This suggests that shocks suppress the short-wavelength features 
by modifying the structure of the aromatic molecules or destroying the smallest grains. Most Seyfert 
nuclei fall on the relationship between aromatic emission and [Ne 11] emission for star-forming galaxies, 
indicating that aromatic-based estimates of the star-formation rate in AGN host galaxies are generally 
reasonable. For the outliers from this relationship, which have small L(7.7 /j,m)/L(11.3 /im) ratios 
and strong H2 emission, the 11.3 /im feature still provides a valid measure of the star- formation rate. 
Subject headings: galaxies: active, galaxies: nuclei, galaxies: Seyfert, galaxies: ISM 



1. INTRODUCTION 

The mid-infrared (mid-IR) aromatic emission features 
are a universal product of star for mation in mo d erate- 
to-high-m etall icity galaxies (e.g ., iRoche et al.l 119911 : 
ILu et al.ll2005lSmith . et al.ll2007ali . Their molecular car- 
riers, often ass umed to be polycyclic aromatic hydro- 
carbons (PAHs.lLeger fe Pugetl ll984; All amandola et al.1 
119851 : iTiel ens 2008) , radiate through IR fluorescence fol- 
lowing v ibrational exc itation by a single ultraviolet (UV) 
photon (Ticlcns 2005) and provide an indirect measure- 
ment of the UV radiation field strength, and therefore the 
star-formation r ate (SFR), that is largely extinction in- 
dependent (e.g..lPeeters et al.ll200llCalzetti et~aLll2007t 
iRieke et a l. 2009). This emission is thought to originate 
in photo-dissociation regions where aromatic molecules 
are s hielded from the ha rsh radiation field near hot stars 
(e.g.. lPovich et alj20 07). These molecules can also be de- 
stroyed by the harder radiatio n field associated with an 
activ e galactic nucleus (AGN, [Voit 1992} iGenzel et al.1 
1998). Nonetheless, aromatic features are readily de- 
tected in many AGNs above IR continua boosted by hot 
dust, and they ha ve been used to probe the SFR in AGN 
host g alaxies (e.g . .lSchweitzer et alj2006t IShi et al.ll2007l 



2009; ILutz et all [20081. Understanding what environ- 
ments destroy or modify these features is important for 
constraining systematic uncertainties in aromatic-based 
estimates of the SFR, and is a key way to probe the 
nature of their molecular carriers, an open issue in our 
u nderstanding of the inters tellar medium. 

iDulev fe~ Williams: ([1981( 1 first suggested that vibra- 
tional modes of aromatic hydrocarbons could produce 
the observed features, which were subsequently identi- 
fied with s pecific C-H and C-C be nding and stretch- 
ing modes (jAllamandola et all [19891. Specifically, the 
6.2 and 7.7 /xm features are produced by C-C stretch- 
ing modes, the 8.6 ^m feature by C-H in-plane bend- 
ing modes, and the 11.3 and 12.7 fim features by C- 
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H out-of-plane bending modes. While these features 
are commonly attributed to PAHs, we use the simpler, 
more general term "aromatic" to avoid making assump- 
tions about the detailed structure of the molecules. It 
is worth noting, for example, that PAH spectra from 
laboratory measurements and quantum chemical calcu- 
lations are unable to match the range of astronomical 
spectra without artificial enhancem ents of the 6 . 2, 7.7 , 
and 8.6 /im feature strengths (e.g.. iLi fc Draind 12001( 1. 
Regardless of this uncertainty associated with uniquely 
matching observed spectra with expectations for specific 
molecules, one can probe the properties of the aromatic 
carriers by measuring the relative strengths of the emis- 
sion features, which are expected to va ry as a function 
of ch arge state (e - g.. IBakes et al.ll200l( l. molecular size 
(e.g.. iDraine fc LiH2001h . and molecular structure (e.g., 
I Vermeii et all 12002( 1 . 

Efforts to study variations in aromatic feature 
strengths outside the Milky Way have focused on star- 
formi ng galaxies (e.g..lSmith et alj|2007at [Galliano et al.l 
[200l iRoseboom et al.l 120091 : lO'Dowd et al.l 12009( 1 . but 
the AGNs included in these studies have shown evidence 
for suppression of shorter wavelength features (e.g., those 
at 6.2, 7.7, and 8.6 /im) relative to th o se at lo nger wave- 
lengths. For example, iSmith et al.l (|2007all studied a 
sample of 59 galaxies fr om the Spitzer Infrar ed Nearby 
Galaxy Survey (SINGS. iKennicutt et al.ll2003l l. of which 
12 have Seyfert nuclei and 20 have low-ionization nuclear 
emission- line regions (LINERs), and they found a ten- 
dency for the systems with reduced 6-8 fj.ui features to be 
associated with low-luminosity AGNs. They speculated 
on possible causes for this behavior, including whether 
the AGN can modify the aromatic grain size distribution 
or serve a s the excitation sourc e for aromatic emission. 
Similarly, lO'Dowd etafl ([2009H studied a sample of 92 
galaxies at z ~ 0.1 from the Spitzer SDSS GALEX Spec- 
troscopic Survey, including eight AGN-dominated and 20 
composite systems, and found that the AGNs exhibited 
lower 7.7/11.3 /jm ratios. They interpreted this behav- 
ior as being consistent with destruction of small aromatic 
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molecules by shocks or X-rays associated with the AGNs, 
but they were unable to distinguish any differences be- 
tween the AGN-dominated and composite objects, nor 
any strong correlation with AGN power. The physical 
slit size at their median redshift is 6 kpc, so there is little 
spatial information. 

In this paper we analyze the aromatic features 
drawing fro m the sample of 89 l o cal S eyfert galaxies 
studied by iDiamon d-Stani c et al.l f|2009T ) . This sam- 
ple is from the revised Sh apley-Ames catalog (RSA, 
ISandage fc Tammannlll987n . and includes every galaxy 
brighter than Bt = 13 that is known to host Seyfert 
activity (jMaiolino fc Riekd 119951: IHo et al.l 11997ft . The 
median distance of the sample is 22 Mpc, so the 3.7" slit 
width of the Sh ort-Low (SL) modu le of the Infrared Spec- 
trograph f IRS. IHouck et aL| | 2004f) on the Spitzer Space 
Telescope (jWerner et al.ll20M ) provides spatial informa- 
tion on scales of a few hundred parsecs, allowing us to iso- 
late nuclear regions distinct from the rest of the galaxy. 
As a result, we are able to probe the effect of AGNs on 
the aromatic features more systematically than has pre- 
viously been done. 

2. DATA 

We gathered data from the Spitzer archive taken with 
the IRS SL module (A = 5.2-14.5 fj,m, R = 64-128) from 
a variety of programs (24, 86, 159, 668, 3247, 3269, 3374, 
30471, 30572, 30577, 30745, 40018, and 50702), as well as 
dedicated data taken for this stu dy (program 40936, PI: 
G. H. Rieke). We use CUBISM (ISmith et alJI2007bD to 
extract one-dimensional spectra from the basic calibrated 
data using 3.6" x 7.2" apertures oriented along the slit 
direction. This aperture size was chosen to isolate the 
nuclear component of the aromatic emission while still 
including a substantial fraction of the diffraction-limited 
beam. We use the calibration for extended sources based 
on the assumption that the regions producing aromatic 
emission are spatially extended, so the extracted spectra 
are in units of surface brightness. We use overlapping 
data in the 7.59-8.42 /im region to scale the SL2/SL3 
orders to the SL1 order when offsets are apparent; these 
offsets are < 10% in all cases. 

We then use a modified version of PAHFIT 
(jSmith et al.l I2007af) to determine the strength of the 
various aromatic features. This spectral fitting pack- 
age includes aromatic features represented by Drude pro- 
files, dust continuum emission represented by modified 
blackbodies at fixed temperatures, fine-structure lines 
and H2 rotational lines represented by Gaussian pro- 
files, starlight represented by T = 5000 K blackbody 
emission, and dust extinction represented by a power- 
law and silicate features. Because Seyfert galaxies ex- 
hibit higher-ionization emission lines, silicate dust emis- 
sion, and hot-dust continuum emission, we additionally 
include a [Ne vi] A7.652 /jm emission line and a silicate 
emission component, both represented by Gaussian pro- 
files, and we use temperatures of 1000, 750, 500, 350, 
225, 150, and 100 K for the dust continuum emission. We 
show example PAHFIT decompositions for three sources 
exhibiting a range in continuum shape and silicate ex- 
tinction in Figure [TJ 
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Fig. 1. — Example PAHFIT spectral decompositions for three 
RSA Seyfert nuclei. The observed spectra are shown in black. The 
blue lines above the continuum correspond to the aromatic fea- 
tures, while red lines correspond to unresolved atomic and molec- 
ular emission lines. The total fit is shown in green, and the dotted 
line indicates the extinction profile. 



3. RESULTS 

The data exhibit a range of signal-to-noise ratios 
(S/N), so we visually inspected all of the nuclear spectra 
to identify those that have clear detections of the relevant 
aromatic features and whose spectra are adequately de- 
scribed by the PAHFIT model. The spectra for these 35 
sources are shown in Figure[2j Since many of the observa- 
tions were executed with the mapping mode of IRS, it is 
also possible to extract spectra for off-nuclear regions in 
some galaxies, allowing for comparison between spectra 
dominated by the active nucleus and spectra dominated 
by H II regions within the same galaxy. We identified 
off-nuclear regions that were covered by the IRS slit and 
had sufficient S/N to detect the relevant aromatic fea- 
tures in 21/35 galaxies. We show a comparison between 
the nuclear and off-nuclear extractions for these galaxies 
in Figure [3l and we compile relevant measurements in 
Tables Q] and E 
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Fig. 2. — Low-resolution 5.2-14.2 (im Spitzer/IRS nuclear spectra for the 35 RSA Seyfert galaxies considered in this study. 
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Fig. [2] — Continued 
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Fig. 3. — Nuclear and off-nuclear spectra for the 21/35 Scyfcrt galaxies where off-nuclear regions were covered by the IRS slit and had 
sufficient S/N to detect the relevant aromatic features. The panel to the right of each spectrum shows the corresponding 3.6" X 7.2" 
extraction region overlaid on the central 1' X 1' of an IRAC 8.0 (im image. 
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Fig. [3] — Continued 
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Fig. 4. — The cumulative distribution of aromatic feature ratios 
for 35 RSA Seyfert nuclei, 21 off-nuclear regions, and 27 SINGS H II 
galaxies. The first three panels show that the 6.2, 7.7, and 8.6 fim 
features are systematically weaker relative to the 11.3 /im feature 
for the Seyfert nuclei than for the off-nuclear regions or the SINGS 
H II galaxies. Panel (c), for example, shows that half of the RSA 
Seyfert nuclei have L(8.6 /im)/L(11.3 /im) ratios < 0.5, whereas 
for the H II galaxies, half have ratios < 0.75. The remaining three 
panels show that the ratios among the 6.2, 7.7, and 8.6 fim features 
show no significant differences between any of the samples. 
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Fig. 5. — The cumulative distribution of the ratio of 6.2, 7.7, and 
8.6 fim features to the 11.3 /an feature for SINGS galaxies with 
Seyfert, LINE R, and H II optical cla ssifications. This illustrates the 
result found bv lSmith et al.l l|2007al ) that the Seyferts and LINERs 
have ratios that are significantly lower than the H II galaxies. The 
apparent difference between SINGS Seyferts and LINERs is not 
statistically significant. 



ISmith et al.1 (|2007al) noted that the LINERs and 
Seyferts in the SINGS sample were offset towards lower 
L(7.7 /tm)/L(11.3 /im) ratios when compared to the H n 
galaxies. We illustrate this result graphically in Figure[5j 
which shows the distribution functions of the ratio of 6.2, 
7.7, and 8.6 /im features to the 11.3 /im feature for galax- 
ies with Seyfert, LINER, and H n optical classifications. 
Both Seyferts (p = 2x 1(T 3 ) and LINERs (p = 0.04) have 
ratios that are significantly lower than the H n galaxies. 
While the SINGS Seyferts have somewhat lower ratios 
than the LINERs, this difference is not statistically sig- 
nificant, and neither sample is statistically distinguish- 
able from the RSA Seyferts. 



3.1. Aromatic Feature Ratio Distributions 

In Figure [H we show the distribution of aromatic fea- 
ture ratios for the 35 RSA Seyfert nuclei and 21 off- 
nuclear regions, as w ell as for 27/59 SINGS galaxies from 
ISmith et alj (J2007a) that have H n nuclear classifications 
(i.e., those that are not Seyferts or LINERs). We find 
that the L(6.2 /im)/L(11.3 /im), L(7.7 /im)/L(11.3 /an), 
and L(8.6 /im)/L(11.3 /im) ratios are systematically 
lower for the Seyfert nuclei than for the off-nuclear re- 
gions or the SINGS H n galaxies. These differences are all 
statistically significant with p < 0.003 based on the two- 
sample Kolmogorov-Smirnov (K-S) test (see Table |3]), 
a non-parametric test that considers the maximum de- 
yiation between two cumulative distrib ution functions 
(jPress et al.lll992t iWall fc Jenk"insl[200l . On the other 
hand, there are no significant differences between these 
feature ratios for the off-nuclear regions and the SINGS 
H II galaxies, so the feature strengths in regions of star 
formation are consistent with being drawn from the same 
parent distribution. Furthermore, the ratios among the 
6.2, 7.7, and 8.6 /im features show no significant differ- 
ences between any of the samples. 



3.2. Trends with Hi emission 

Inspection of Figure [5] reveals that several Seyferts 
with small L(7.7 /im)/L(11.3 fim) aromatic feature ratios 
also exhibit a strong H2 S(3) rotational line at 9.67 /im 
(e.g., NGC4501, NGC5194.) To investigate this behav- 
ior, we plot the strength of the H2 S(3) line, normal- 
ized to the strength of the aromatic features, as a func- 
tion of the L(7.7 /im)/L(11.3 /im) ratio in Figure |5J 
We find a strong anti-correlation in this plot such that 
sources with the smallest L(7.7 /tm)/L(11.3 /im) ratios 
also have the strongest H2 emission. The Spearman's 
p rank correlation coefficient is —0.78 with a probabil- 
ity p = 3 x 10~ 8 of no correlation, while Kendall's r 
is —0.62 with p = 1 x 10~ 7 ; these non-parametric tests 
consider the agreement between the ranks of quantities in 
pairs of measurements ([Press et al.l ll992; Wall fc Jenkins! 
2003), with coefficient values ranging fro m -1 (perfect 
disagr eement) to 1 (perfect agreement). iRoussel et al.l 
(2007) found that H2 rotational lines scale tightly with 
the aromatic features for SINGS H 11 galaxies, but that 
Seyferts and LINERs often exhibit excess H2 emission, 
which they attribute to shocks. We explore the hypoth- 
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Fig. 6. — The relationship between the strength of the H2 S(3) 
rotational line, normalized to the strength of the aromatic fea- 
tures, and the L(7.7 /im)/L(11.3 fim) ratio for RSA Seyfert nu- 
clei. The sources with small L(7.7 /tm)/L(11.3 /im) ratios also 
exhibit strong H2 emission. The most extreme sources with 
L(7.7 /xm)/L(11.3 fim) < 1.6 are highlighted with filled circles. 



esis that shocks cause both the excess H2 emission and 
the anomalous aromatic ratios for AGNs in Section l4~4l 
Among the sources excluded from the above analysis 
due to a lack of 6.2, 7.7, or 8.6 fim aromatic feature de- 
tections, there are a significant number with clearly de- 
tected 11.3 /im features and H 2 S(3) lines. In Figure El 
we show the nuclear spectra for a dozen of these sources, 
sorted by the equivalent width of the 11.3 /im feature. 
These spectra exhibit the small L(7.7 ^m)/L(11.3 fim) 
ratios and strong H 2 S(3) lines characteristic of sources in 
the top-left of Figure El Due to uncertainties associated 
with estimating the strength of weak, broad features and 
determining robust upper limits (e.g., proper continuum 
placement), we do not include any of these sources in 
our subsequent analysis. However, their behavior is con- 
sistent with that in Figure |5] and supports the reality of 
the trend between aromatic feature characteristics and 
H 2 line strength. 

3.3. Evidence for extinction of aromatic features 

The sources with the largest L(7.7 /im)/L(11.3 fim) ra- 
tios, NGC4945 and NGC3079, also have the strongest sil- 
icate absorption features. This suggests that the 11.3 fim 
feature is being significantly attenuated, consistent with 
previous res ults for starburst and luminous infrared 
galax ies (e.g.. lBrandl et al] |2006; Pcrcira-Santaell a et al.l 
I2010D . and implies that a significant fraction of the 
silicate-absorbing material is extended relative to the re- 
gions that produce the aromatic features. Although the 
aromatic feature measurements in PAHFIT are corrected 
for extinction, in cases as extreme as these two galaxies 
the resulting feature strengths are highly uncertain. For 
all other galaxies in our sample, the inferred extinctions 
are < 50% for all features. 
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Fig. 7. — Nuclear spectra for 12 additional RSA Seyfert nuclei 
that exhibit small L(7.7 /^m)/L(11.3 fim) ratios and strong H2 S(3) 
lines, but were excluded from the sample due to a lack of 6.2, 7.7, or 
8.6 fim aromatic feature detections. The spectra are sorted from 
top to bottom by the equivalent width of the 11.3 fim aromatic 
feature. The wavelengths of the 9.67 fim H2 S(3) line and the 
11.3 fim aromatic feature are marked by dotted lines. 



4. DISCUSSION 

The result that Seyfert galaxies exhibit weak 6.2, 7.7, 
and 8.6 fim aromatic features relative to the 11.3 fim 
feature could be explained by radiative or mechanical 
processing of the molecular carriers by the active nucleus. 
Here we explore the relevant physical and chemical effects 
that could modify the observed feature strengths. 

4.1. Ionization Balance 



Previous experimental (e.g., [Szczepanski fc Valalll993l : 
Hud gins &: Allamandolal Il995h an d theoretical (e.g. 
331: iLane 



DeFrees et al.l 119931 : lLanghofi1H996l ) work on PAHs has 



shown that the C-C stretching modes that produce the 
6.2 and 7.7 fim features, as well as the C-H in-plane 
bending modes that produce the 8.6 fim feature, are 
more efficiently excited in ionized molecules. The ra- 
tios of these features to the 11.3 fim feature, which is 
produced by C-H out-of-plane bendi ng modes, are lower 
for n eutral molecules (see Figure 1 of lAllamandola et al.l 
1999). The fraction of ionized aromatic molecules is set 
by the balance between ionization and recombination, 
which depends on the UV radiation field density (Go), 
the gas temperature (T), and the electron density (n e ) 
a ccording to G n T 1/2 /n e (jBakes fc Tielenslfl99l . 

iGalliano et al.l (|2008l ) argued that the variations in aro- 
matic feature ratios for a heterogeneous sample of 50 
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objects (including Galactic regions, Magellanic H n re- 
gions, and galaxies, as well as spatially resolved regions 
within seven of those objects) are controlled by this ion- 
ization balance. Similar to our results, they found that 
the relative strengths of 6.2, 7.7, and 8.6 /jm features 
showed little variation, while the ratios between these 
features and the 11.3 /im feature varied by an order of 
magnitude. This hypothesis is also supp orted by ob- 
servat ions of Galactic reflection neb ulae by Uoblin et al.l 
(1996) and lBregman fc Temil (|2005l ). who found decreas- 
ing L(8.6 /tm)/L(11.3 /im) and L(7.7 /rai)/L(11.3 fxm) 
ratios as a function of distance from the ionizing source, 
consistent with an increasing neutral fraction. 

To compare with model expectations for ion- 
ized and neutral aromatic molecules, we plot 
L(11.3 /im)/L(7.7 /on) v. L(6.2 /tm)/L(7.7 /on) 
for the Seyfert nuclei, off-nuclear sources, and SINGS 
H II galaxies in Figure [HI This c an be compared 
with Figure 1 6 of | Draine fell (l200l and Fi gure 5 of 
IQ'Dowd et all (120091 7 although we use a condensed plot 
range. We find that a number of Seyferts lie beyond the 
range of model predictions, even for completely neutral 
aromatic molecules; these are the 6/35 Seyferts with 
L(11.3 /zm)/L(7.7 //m) > 0.6: NGC5194, NGC4501, 
NGC4639, NGC1433, NGC2639, and NGC5005. While 
such L(11.3 /im)/L(7.7 /im) ratios could be produced 
by large (> 200 C atoms) neutral molecules, they would 
be expected to have L(6.2 /tm)/L(7.7 fim) < 0.25, 
which is inconsistent with the data. Sim ilar extreme 
aroma tic band strengths were observed by iReach et al.1 
(120001) for the quies cent molecular cloud SMC Bl No. 
1, and iLi fc Draind {2002) were unable to reproduce 
the observed band ratios even with completely neutral 
grains. 

The above comparison is for a single Milky Way- 
based model, and laboratory studies have found larger 
L(11.3 /tm)/L(7.7 /tm) ratios for neutral PAHs, but it 
does illustrate the difficulty in explaining our results for 
Seyfert galaxies in terms of a low ionized fraction. Fur- 
thermore, under the assumption that aromatic features 
are produced by star formation (see Section I4.5j) . the 
temperatures and densities of the aromatic-emitting re- 
gions should be typical of PDRs, whereas the UV radi- 
ation field would likely be enhanced by the AGN. This 
implies that the ionized fraction would be higher, not 
lower. Thus, ionization balance arguments appear un- 
able to explain the behavior of the aromatic features 
around AGNs. 

4.2. Grain Size 

Smaller aromatic molecules contribute preferentially 
to th e shorter- wavelength features (e.g., iSchutte et al.1 
1993), but they are subject to photodestruction by the 
UV radiation field and collis ional destruction by shocks. 
Based on laboratory studies. IJochims et al.l (|1994f ) found 
a critical size of 30-40 C atoms, belo w which PAH s would 
mainly be photodissociated, while lAllain et al.l (1996) 
suggested a larg er critical value of 5 C atoms based 
on their models. iLe Page et al.l (|2003l ) agreed that small 
PAHs with 15-20 C atoms or fewer would be destroyed 
in most environments, but their models indicated that 
PAHs in the 20-30 C atom range may survive, albeit 
with most of their peripheral H atoms stripped away, 
while larger PAHs would survive with their H atoms in- 
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Fig. 8. — The relative strengths of the 6.2, 7.7, and 11.3 fim fea- 
tures for RSA Seyfert nuclei, off-nuclear regio ns, and SINGS H I I 
galaxies compared with model predictions from[Drainc & Lj| 1120011 ) 
for neutral and ionized PAHs. The dashed lines correspond to pre- 
dictions for completely neutral and completely ionized molecules; 
the permitted region of the diagram is bounded by these two lines. 
The arrows illustrate the effects of increasing grain size and in- 
creasing ionization on the aromatic feature ratios. The Seyferts 
highlighted as filled circles in Figure [6] all lie beyond the range of 
model predictions, even for completely neutral molecules. 



tact. iMicelotta et al.l (|2010al) found that PAHs with 50 C 
atoms would not survive in shocks with velocities greater 
than 100 km s _1 , while PAHs with 200 C atoms would 
be destroyed by shocks with velocities above 125 km s _1 . 

Destruction of the smallest molecules is expected to 
result in the 6.2 and 7.7 /im features being suppressed 
relative to the 11.3 /im feature, as well as the 6.2 /tm 
featu re being suppressed relative to the 7.7 /tm feature 
fe.g.. IDraine fcIJl200l iGalliano et alJ 12008ft . The for- 
mer effect is clearly seen in Figure 21 but the latter is not. 
Thus, the hypothesis that small-grain destruction can ex- 
plain the observed ratios is only tenable if the molecules 
that produce the 6.2, 7.7, and 8.6 /im features are de- 
stroyed with similar efficiency, which is inconsistent with 
existing models. 

4.3. Hydrogenation and Molecular Structure 

The level of hydrogenation of the aromatic molecules 
will affect the number of C-H bonds and therefore 
the relative strength of the C-H and C-C vibra- 
tional modes. An increas e in th e C-H/C-C ratio was 
proposed by IReach et al.l (j2000P i to explain the large 
L(11.3 /zm) /L (7. 7 urn) ratio o bserv ed in SMC Bl No 
1, although iDrainefe Lil (|2001[ ) and ILi fc Draind (|2002f ) 
point out that PAHs with > 30 C atoms are already 
expected to be fully hydrogenated. Some range in C- 
H/C-C ratios, even for fully hydrogenated molecules, 
is facilitated by the structure of the C skeleton, which 
can be compact with more C-C bonds or open with 
more C-H bonds (e.g., perico ndensed PAHs v. cata- 
condensed PAHs, iTielensi [20051 ). The structure also af- 
fects the number of adjacent C-H groups per aromatic 
ring, and therefore the relative strengths of the 11.3 /tm 
feature, which is produced by solo C-H bonds, and the 
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Fig. 9. — The cumulative distribution of L(12.7 /im)/L(11.3 £tm) 
ratios for RSA Seyfert nuclei, off- nuclear regions, and SINGS H II 
galaxies. The result that Seyfert nuclei exhibit significantly smaller 
ratios suggests aromatic molecules that have fewer adjacent C— H 
groups. 



12.7 / xm feature, which is produced by C-H multiplcts 
(e.g., iHonv et al.l [2001). For example, based on the 
large L(12.7 /mi)/L (11.3 /mi) ratio for SMC Bl No. 1, 
iVermeii et al.l (2002) argued for a compact structure with 
a higher incidence of C-H multiplets. 

To investigate such behavior, we plot the 
L(12.7 /mi)/L(11.3 /mi) ratios for Seyfert nuclei, 
off- nuclear regions, and SINGS H II galaxies in Figure [SI 
The Seyfert nuclei exhibit significantly smaller ratios 
(p < 0.001), while the ratios for off-nuclear regions 
and SINGS H n galaxies are not distinguishable. This 
implies that the aromatic molecules in Seyfert nuclei 
may have fewer C-H multiplets. Thus a scenario where 
AGN processing or environment results in open, uneven 
molecular structures with higher C-H/C-C ratios and 
fewer adjacent C-H groups could qualitatively explain 
the observed 6-13 /im aromatic spectra. 
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Fig. 10. — The relationship between the aromatic features and 
the [Ne II] and [O iv] emission lines. The filled circles correspond 
to the RSA Seyferts defined in Figure [6] that have the smallest 
L(7.7 /im)/L(11.3 /im) ratios. The strong correlation with [Ne III, 
which traces star formation, and the weak correlation with [O iv], 
which traces AGN activity, implies that the aromatic features are 
primarily associated with star formation. Most of the Seyfert nuclei 
lie on the relationship between aromatic feature and [Ne n] emission 
for H II galaxies; the only outliers are among the sources highlighted 
with filled circles, which have extreme aromatic feature ratios (see 
Figures \E\ and [8} . The aromatic feature and [Ne II] emission values 
are in surface bri ghtness units (W m~ 2 sr —1 ) , while the [O iv] 
values, taken from Diamond-Stanic et alJ (2009), are in flux units 
(W m" 2 ). 



4.4. The role of AGN -driven shocks 

As presented in Section 13.21 and Figure El the Seyfert 
galaxies with the smallest L(7.7 /mi)/L(11.3 /urn) aro- 
matic feature ratios also exhibit the strongest H2 S(3) 
emission, which probes hot molecular gas (upper level 
temperature 2500 K). The incidence of this excess H2 
emission does not scale with AGN luminosity, indi- 
cating that shock e xcitation is more im portant than 
X-ray heating (e.g., iRoussel et al.l 120071 ). A connec- 
tion between shock-heated, ^-emitting gas and small 
L(7.7 /mi)/L(11.3 /mi) ratios was found by lOgle et al. 
( 20071) for the radio galaxy 3C 326 and bv lGuillard et ah 
(2010) for Stephan's Quintet, a compact group of in- 
teracting galax ies exh ibiting a larg e-scale shock (e.g., 
Appleton et"aLl 120061 : iCluver et alT l2010f ). Similarly, 
Kan eda et al. 120081) " found strong H2 emission and small 
L(7.7 /tm)/L(11.3 /im) ratios in a sample of local el- 
liptical galaxies, man y of which host low-luminosity 
AGNs. More recently, iVega et al.1 (|2010f ) affirmed this 



result for a sample of four early-type galaxies classi- 
fied as LINERs, and they argued that shock process- 
ing of aromatic molecules may be responsible for the 
observed behavior. As discussed above, the Seyferts 
and LINERs in the SINGS sa mple also exhibit sm aller 
L(7.7 /mi)/L(11.3 urn ) ratios (iSmith et al.1 12007a!) and 
stronger H2 emission (jRoussel et al.1 l2007f ) than do the 
H 11 galaxies. 

Shocks are expected to have profound impacts on in- 
terstellar dust via shattering in grain-g rain collisions and 
sputt ering in ion-grain collisions fe.g.. Uones et al.lll994 
1996). Aromatic features are nonetheless observed in 
the shocked environments of supern ova remnants (e.g., 
iTappe e t~a l. 2006; Reach ct al. 2006) and galactic winds 
(e.g., iTacconi-Garman et al.l 12005c lEngelbracht et all 
2006). The observed emission may come from en- 
trained clumps that are not fully exposed to the sho ck 
or the hot, po s t-shock gas ([Micelotta et al.1 l2010allbl ) . 
iMicelotta et al.1 ([2010a) study the processing of small 
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Fig. 11. — The relationship between aromatic feature and [Ne n] 
emission as a function of the L(7.7 /im)/L(11.3 fim) ratio. The 
filled circles correspond to the RSA Seyferts defined in Figure [6] 
that have the smallest L(7.7 /im)/L(11.3 fim) ratios. This con- 
firms that the sources with suppressed aromatic features, relative 
to [Ne n], have the smallest L(7.7 ^m)/L(11.3 ^tm) ratios. 



carbon grains (Nc < 200, corresponding to aromatic 
molecules) by interstellar shocks and find that their 
molecular structure is severely denatured for shock ve- 
locities of 75-100 km s _1 and they are completely de- 
stroyed when v > 125 km s" 1 . The effect of this shock 
processing on the observed aromatic feature ratios is not 
known. A possibility that could explain the association of 
modified aromatic feature ratios with strong H2 emission 
is that shocks may leave open, uneven structures in the 
surviving aromatic molecules. We note that AGN-driven 
shocks, if responsible for the observed behavior, do not 
strongly suppress the 11.3 pm feature (see Section 14751) . 

4.5. Could the aromatic features be excited by the 
AGN? 

iSmith et all (|2007aD speculated that the AGNs could 
directly excite aromatic emission. If this were the case, 
SFRs estimated from aromatic features would be over- 
estimated due to this AGN contribution. To investi- 
gate the relationship between star formation rate, AGN 
luminosity, and aromatic feature strength, we plot the 
fluxes of the [Ne II] and [O iv] emission lines versus 
those of the 7.7 (im and 11.3 pm aromatic features in 
Figure [TO] The [Ne 11] line has an ionization potential 
of 21 eV and is a reasonable tracer of the SFR (e.g., 
IHo fc Ketoll2007l ). while the [O iv] line has an ionization 
potential of 55 eV and traces the AGN intrinsic lumi- 
nosity (e.g..lMelendez et all 2008; Di amond-Stanic et all 
2009; Rigby et al.l l2009ri . Figure [TU] shows the strong cor- 
relation between [Ne n] and aromatic feature strength 
for RSA Seyferts (Spearman's p = 0.93), which matches 
the relationship for SINGS H II galaxies, and it shows 
the weak correspondence between [O iv] and aromatic 
feature strength (Spearman's p = 0.39). This confirms 
that the aromatic features are primarily tracing star- 
formation activity. 

The Seyferts that are outliers in the [Ne n] -aromatic 
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Fig. 12. — The relationship between [Ne II] emission and the 
7.7 and 11.3 fim aromatic features. The filled circles correspond 
to the RSA Seyferts defined in Figure [6] that have the smallest 
L(7.7 /im)/L(11.3 /an) ratios. While the 7.7 fim feature can be 
strongly suppressed, the 11.3 /im feature is still a robust tracer 
of the SFR. The solid line in the bottom panel corresponds to the 
median ratio Lj]\j n ]/Ln 3 = 0.12, and the dotted lines correspond 
to factors of two above and below this median value. Scatter in this 
ratio is expected because [Ne 11] traces somewhat younger stellar 
populations than do the aromatic features. All values are in surface 
brightness units (W m — 2 sr _1 ). 



feature relationship have weak aromatic features, and we 
show in Figure [TT] that these correspond to the sources 
with the smallest L(7.7 /im)/L(11.3 pm) ratios. There 
are no examples with stronger aromatic features as might 
be expected if the AGN were exciting additional emis- 
sion. The three obvious outliers, NGC2639, NGC4501, 
and NGC5194, all have [O iv]/[Ne n] < 0.25, imply- 
ing t hat the AGN con tribution to [Ne n] is < 10% 
(e.g., iSturm et all I2002T) . We note that the incidence 
of modified aromatic spectra does not show a depen- 
dence on AGN lu minosity, confirming the results of 
IBaum et~a l (2010), who found no correlation between 
the L(6.2 /im)/L(11.3 pm) ratio and [Ne v] luminosity. 

4.6. Use of Aromatic Features to Determine SFRs 

Several studi es (e.g . JSchweitzer et a l. 2006; Lutz et all 
120081: IShi et"alll2009l ) have used the 6.2 and 7.7 pm aro- 
matic features to measure the SFRs in AGN host galax- 
ies. The result that some AGNs exhibit suppressed short- 
wavelength aromatic features (e.g., the outliers in Fig- 
ure llip suggests that such SFR measurements may be un- 
derestimated. To determine whether the 11.3 pm feature 
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is robust to such effects, we plot separately the relation- 
ships between [Ne n] and the 7.7 and 11.3 /mi features in 
Figure QU We find that almost all of the RSA Seyferts, 
including those with anomalously high L[Nenl/L7.7 val- 
ues, are within a factor of two of the median value 
L[Nen]/Lii.3 = 0.12. Scatter in this ratio is expected 
as a function of the age of the stellar population because 
21 eV photons from young stars (< 10 Myr) are required 
to produce [Ne n], while somewhat older stars can pro- 
duce 6-13. 6 eV UV photons t h at excite aromatic emis- 
sion (e.g.. iPeeters et al.l l2004t iDiaz-Santos et all 120101 : 
iPereira-Santaella et al.l 12010( 1. Silicate absorption will 
tend to increase the observed ratio, but this is only a sig- 
nificant effect for sources like NGC4945 and NGC3079 
(see Section !3-3[) . While SFR estimates based on the 
11.3 /im feature are still subject to the uncertaintie s 
that apply to H n galaxies (e.g., iSmith et al.l l2007al ). 
such measurements for AGN hosts appear to be robust 
to the effects of AGN- and shock-processing of aromatic 
molecules. 

5. CONCLUSIONS 

We have shown that the relative strengths of the mid- 
IR aromatic features for Seyfert galaxies differ signifi- 
cantly from those for star-forming galaxies, with the 6.2, 
7.7, and 8.6 /im features being suppressed relative to 
the 11.3 /an feature in Seyferts. The sources with the 
smallest L(7.7 /mi)/L(11.3 /im) aromatic feature ratios 
also exhibit the strongest H2 S(3) rotational lines, which 
likely trace shocked gas (see Figure [6|). We explore the 
relevant physical and chemical effects that could pro- 
duce the observed aromatic spectra. An enhanced frac- 
tion of neutral aromatic molecules could produce quali- 
tatively similar behavior, but the observed ratios lie be- 



yond model predictions for completely neutral molecules 
and the presence of an AGN would be expected to in- 
crease the level of ionization rather than reduce it. De- 
struction of the smallest aromatic molecules could ex- 
plain the suppression of shorter wavelength features, but 
the expected variations in the relative strengths of the 
6.2, 7.7, and 8.6 /im features are not seen. A mod- 
ification of the molecular structure that enhances the 
C-H/C-C ratio could reproduce the observed behavior, 
and an open C skeleton with fewer adjacent C-H groups 
would furthermore explain the reduced strength of the 
12.7 /im feature. Given the connection between strong 
H2 emission and modified aromatic ratios, we speculate 
that shock processing could produce such structures. Fi- 
nally, we show that the aromatic features correlate well 
with [Ne 11] (i.e., star formation) but not with [O iv] 
(i.e., AGN luminosity), indicating that AGN excitation 
of aromatic emission is not significant and that aromatic- 
based estimates of the SFR are generally reasonable. 
There are a few outliers with strong H2 emission, small 
L(7.7 /mi)/L(11.3 /mi) ratios, and small aromatic/[Ne 11] 
ratios, but for these sources the 11.3 /mi feature is still a 
reasonably robust tracer of the SFR. 
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TABLE 1 

Nuclear Measurements 
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Note. — Measurements arc in units of W m 2 sr 
a Consists of sub-features at 7.42, 7.60, and 7.85 (xm. 
b Consists of sub-features at 11.23 and 11.33 fim. 
c Consists of sub-features at 12.62 and 12.69 fim. 



The Effect of AGNs on Aromatic Features 



TABLE 2 

Off-Nuclear Measurements 
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32e-07 


4 


.59+0 


,25e-07 
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13e-06 


8. 


.07+0 


25e-07 
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Note. — Measurements arc in units of W m 2 sr 
a Consists of sub-features at 7.42, 7.60, and 7.85 /im. 
b Consists of sub-features at 11.23 and 11.33 /^m. 
c Consists of sub-features at 12.62 and 12.69 /xm. 



TABLE 3 
Statistical Tests 



ratio Seyferts v. SINGS Seyferts v. off nuclear SINGS v. off nuclear 



6/11 5 x 10~ 4 5 x 10 -5 0.682 

7/11 0.003 0.001 0.074 

8/11 9 x 10~ 4 2 x 10~ 4 0.063 

6/7 0.230 0.447 0.888 

6/8 0.489 0.347 0.689 

7/8 0.108 0.303 0.374 



Note. — 



Values correspond to probabilities from two-sample K-S tests. 



